AD-A058  050 


UNCLASSIFIED 


KEARNEY  (A  T)  INC  CHICAGO  IL  CAYWOOO-SCHILLER  DIV  F/G 

DYNAMIC  SAM  ENDGAME  MOOEL  ANALYST’S  MANUAL.  VOLUME  I.(U) 

SEP  77  M A DLOOGATCH.  R H ROSE.  0 S KLUK  N62269-76-C-0386 

CSD-77-200-VOL-1  NL 


F/G  16/4.2 


1 <*  1 

ADA 

058050 

■ 

M 

g 

P 

N 

S 

END 

DATE 

FILMED 

10  78 

DOC 

ADA058050 


'H 


\ [ 

‘ — * ;J  L . 


/ 


DYNAMIC  SAM  ENDGAME  MODEL 
ANALYST'S  MANUAL 
N62269-76-C-0386 

SEPTEMBER  1977 


Kearney  Management  Consultants 


CAYWOOD- SCHILLER  DIVISION 
A.  T.  KEARNEY,  INC. 

100  South  Wacker  Drive 
Chicago,  Illinois  60606 


t DYNAMIC  SAM  ENDGAME  MODEL 

ANALYST’S  MANUAL 

> 

> FINAL  REPORT 

Under 

Contract  N62269-76-C-0386 


[ 


Prepared  for 

NAVAL  AIR  DEVELOPMENT  CENTER 
WARMINSTER,  PENNSYLVANIA  18974 

78  0 8 25  041 


TABLE  OF  CONTENTS 


Section 

1.  Introduction  

2.  Detailed  Description  of 
Endgame  Additions  to  DSAMAM 

2 . 1 Logical  Flow  

2.2  Coordinate  Systems 

and  Transformations  . . . 

2.3  Fuzing  

2.3.1  Conical  Radar  Fuze 

2.3.2  Contact  Fuze  . . . 

2.4  Single-Shot  Component  Kill 

2.4.1  Direct-Hit  Kill  . 

2.4.2  Blast  Kill  .... 

2.4.3  Fragment  Kill  . . 

2.5  Overall  Aircraft  Kill  . . 


Page 

1 

2 

2 

6 

9 

9 

12 

13 

14 
14 
14 
22 


□ □ 


LIST  OF  FIGURES /TABLES 


Figure  No. 

2-1  Logic  for  PK  Determination  

2-2  Static  and  Dynamic  Fragment  Emission 

Table  No. 


2.1 


Criteria  for  Aspects  Struck 


FOREWORD 


The  activities  described  in  this  report  in  modifying  the 
Dynamic  SAM  Model  to  include  an  endgame  capability  were 
performed  during  the  period  July  1976  to  September  1977 
under  Contract  N62269-76-C-0386  for  the  Naval  Air  Development 
Center,  Warminster,  Pennsylvania.  The  purpose  of  this  effort 
was  to  develop  a means  for  evaluating  kill  probabilities 
associated  with  trajectories  and  missile/aircraft  intercept 
conditions  generated  by  the  NADC  Dynamic  SAM  Model.  This  is 
) the  first  of  two  volumes  constituting  the  final  report.  The 

work  was  performed  by  R.  H.  Rose,  M.  A.  Dloogatch,  and 
D.  S.  Kluk  of  the  Caywood- Schiller  Division  of  A.  T.  Kearney. 
t 

This  volume  includes  a summary  of  the  work,  and  an  analytic 
description  of  the  model. 

» 
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1 . INTRODUCTION 
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The  Dynamic  SAM-Aircraft  Model  (DSAMAM)  was  originally 
designed  and  implemented  by  Autonetics  in  1969.  The  intent 
was  to  provide  a tool  for  evaluating  the  effectiveness  of  a 
maneuvering  aircraft  versus  a command  guided  surface-to-air 
missile.  The  output  of  this  model  consisted  of  the  inter- 
cept geometry  for  each  missile.  No  attempt  was  made  to 
determine  a fuzing  or  detonation  point  or  a kill  probability. 

This  original  model  was  programmed  in  FORTRAN  IV  to  run  on  an 
IBM  360  computer.  At  some  later  date  the  necessary  changes 
were  made  to  allow  the  model  to  be  run  on  the  CDC  computer 
at  NADC.  At  this  time  a modification  was  made  so  that  a kill 
probability  was  calculated  based  on  the  assumption  of  detona- 
tion at  the  point  of  closest  approach.  This  PK  was  based  on 
aircraft  lethal  radius  inputs  and  a specified  CEP  for  the 
missile . 

The  objective  of  the  current  study  was  to  provide  DSAMAM  the 
capability  of  simulating  the  fuzing  characteristics  and  ter- 
minal effects  of  both  the  SA-2  and  SA-3  proximity-fuzed 
warheads . 

Several  modular  subroutines  were  built  into  DSAMAM  in  order 
to  complete  this  objective.  These  allow  for  a conical-pattern 
radar  fuze,  a backup  contact  fuze,  calculation  of  blast  and 
fragment  kill  for  each  component  of  the  aircraft,  and  a cumu- 
lative kill  probability  for  the  entire  aircraft. 
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2.  DETAILED  DESCRIPTION  OF  ENDGAME 
ADDITIONS  TO  DSAMAM 

The  following  description  assumes  that  the  reader  is  familiar 

with  Volumes  I-IV  of  the  Final  Report  on  the  Dynamic  SAM- 

Aircraft  Study,  Report  No.  C9-2336/120,  Autonetics,  North 

American  Rockwell  Corporation. 


2 . 1 Logical  Flow 

The  flow  chart  of  Figure  2-1  depicts  the  logic  involved  in 
determining  whether  or  not  fuzing  has  occurred  at  each  time 
step,  and  if  it  has,  calculating  kill  probabilities  for  each 
component.  The  notation  used  in  these  flow  charts  is  as 
follows : 


DHK 

DHNK 

Pd(CF) 

pd(RF) 

Tb(CF) 

Tb(RF) 

T (DHK) 

T (RF) 

T (DHNK) 
Pk(T) 

CONFUZE 

SPRAY 

STAB 


direct  hit  with  kill 

direct  hit  without  kill 

probability  that  contact  fuze  is  dud 

probability  that  radar  fuze  is  dud 

time  of  detonation  due  to  contact  fuze 

time  of  detonation  due  to  radar  fuze 

time  at  which  direct  hit  with  kill  occurs 

time  at  which  radar  fuzing  occurs 

time  at  which  direct  hit  without  kill  occurs 

probability  of  component  kill  due  to  fragment 
spray  given  detonation  at  time  T 

subroutine  which  determines  whether  or  not  the 
conical  radar  fuze  has  been  activated 

subroutine  which  calculates  the  component  kill 
probabilities  resulting  from  fragment  damage 
due  to  an  explosion  at  some  specified  time 

subroutine  which  determines  whether  any  of  a set 
of  direct-hit  ellipsoids  has  been  pierced  by  the 
path  of  the  missile  during  the  current  time  step 
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FIGURE  2-1 

LOGIC  FOR  PK  DETERMINATION 


DAK  n,-ZINC\JL 
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CALL  SPRAY 

CALL  SPRAY 

(T-Tb(CD) 

(T*Tb(CF)) 

SET  COMPONENT 
PK ' S TO 

(l-Pn(CF))*F*(T.(CF)) 


SET  COMPONENT 
PK’S  TO  1 


SET  COMPONENT 
PK'S  TO 

(1-PD(CF))*FK(TE(CF))+PD(CF) 


COMBINE  COMPONENT 
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A/C  KILL  PROBABILITY 
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FIGURE  2-1  (continued) 
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FIGURE  2-1  (continued) 
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The  logic  takes  into  account  the  possibility  of  the  missile 
possessing  either  a radar  fuze  or  a contact  fuze  or  both. 

Each  fuze  is  considered  to  have  a dud  probability  associated 
with  it.  This  dud  probability  would  be  supplied  as  an  input 
by  the  user.  Thus,  for  the  case  of  a missile  which  has  no 
contact  fuze,  the  contact  fuze  dud  probability  would  simply 
be  set  to  1.0. 

An  assumption  which  is  implicit  in  the  logic  used  is  that  the 
radar  fuze  will  not  be  activated  after  a direct  hit  has  occurred. 

2 . 2 Coordinate  Systems  and  Transformations 
In  order  to  avoid  any  confusion  arising  from  differing  sign 
conventions,  the  coordinate  systems  used  in  the  derivations 
of  the  endgame  analysis  were  chosen  to  be  identical  to  those 
defined  in  Autonetics ' Final  Report  on  the  Dynamic  SAM-Aircraf t 
Study.  The  three  rectangular  coordinate  systems  employed  are: 

1)  An  inertial  coordinate  system  (defined  as  the  "aircraft 
state  vector  coordinate  system,"  page  2-10,  Volume  II;  see 
Figure  4-1,  page  4-7,  Volume  III).  This  coordinate  system  has 
its  origin  at  the  SAM  site.  The  Z-axis  is  perpendicular  to  the 
surface  of  the  earth  at  the  SAM  site  and  positive  upwards.  The 
X-axis  and  the  Y-axis  lie  in  the  plane  perpendicular  to  the  Z- 
axis  and  are  defined  so  as  to  create  a right-handed  system. 

2)  An  aircraft  body  coordinate  system  (defined,  page  3-28, 
Volume  II).  The  X-axis  of  this  coordinate  system  is  parallel 

to  the  thrust  vector  of  the  aircraft.  The  Y-axis  is  perpendicular 
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to  the  X-axis  and  positive  in  the  direction  of  the  left  wing. 
The  Z-axis  is  positive  upward  through  the  pilot's  head. 

Because  in  the  basic  DSAMAM  program  the  aircraft  is  essentially 
considered  to  be  a single  point,  the  location  of  the  origin 
relative  to  the  components  of  the  aircraft  is  not  specified  in 
the  definition  of  the  coordinate  system.  However,  since  the 
endgame  addition  represents  the  aircraft  in  a variety  of 
ways  (e.g.,  a set  of  points  or  a set  of  ellipsoids)  the  loca- 
tion of  the  origin  must  be  specified  for  the  aircraft  body 
coordinate  system.  The  choice  is  arbitrary;  however,  the 
user  will  probably  wish  the  point  to  correspond  to  some  natural 
feature  of  the  aircraft  (e.g.,  nose,  center  of  gravity,  pilot). 

3)  A missile  body  coordinate  system  (defined,  page  2-19, 
Volume  III;  see  Figure  2-6,  Volume  III).  The  X-axis  in  the 
missile  body  syrtem  corresponds  to  the  longitudinal  axis  of 
the  missile  and  is  positive  towards  the  nose.  The  Y-axis  is 
perpendicular  to  the  X-axis  and  lies  within  a plane  parallel  to 
the  inertial  XY-plane.  The  Y-axis  is  positive  to  the  right 
of  an  observer  sitting  on  the  missile  and  facing  forward.  The 
Z-axis  is  perpendicular  to  the  other  two  and  chosen  so  as  to 
form  a right-handed  coordinate  system.  Thus,  for  a missile 
flying  level,  the  Z-axis  is  positive  downwards. 

In  order  to  facilitate  conversion  between  these  coordinate 
systems,  three  transformation  matrices  are  computed  and  stored 
in  COMMON  storage.  This  is  done  in  subroutine  MISSIL.  These 
matrices  are  named  TMI , TIA,  and  TAM. 
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The  first  matrix,  TMI,  is  so  named  because  multiplication  of 
a column  vector  by  this  matrix  will  rotate  the  vector  from  a 
coordinate  system  based  on  the  Missile  body  axes  to  the 
Inertial  coordinate  system.  Similarly,  multiplication  by 
TIA  rotates  a column  vector  between  the  Inertial  coordinate 
system  and  a system  based  on  the  Aircraft  stability  axes, 
and  multiplication  by  TAM  rotates  a column  vector  from  the 
Aircraft  system  to  the  Missile  system. 

These  matrices  are  calculated  as  follows: 


/ cos(eM)  cosO'M)  -sin(GM)  cos(eM)  sinC1^) 

TMI  = j-sin(eM)  cos(TM)  -cos(eM)  -sin(6M)  sin^) 


sin(4-M) 


■cos^) 


where  9^  and  4^  are  the  elevation  and  azimuth,  respectively, 
of  the  missile  body  axis  in  the  inertial  coordinate  system 
(defined,  page  2-19,  Volume  III  of  the  Final  Report).  4>M 
is  positive  when  the  missile  is  above  the  inertial  XY  plane 
and  9^  is  positive  when  the  positive  X-axis  in  the  missile 
body  system  is  to  the  right  of  the  positive  X-axis  in  the 
inertial  system. 


A 


where  T,  aIf  and  dj  are  the  unit  vectors  along  the  aircraft 
body  axes  in  the  inertial  coordinate  system  (defined,  page 
3-22,  Volume  II). 
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TAM  - (TMI)T  (TIA)T  , 

where  the  superscript  T signifies  the  transpose  of  the  matrix. 

2 . 3 Fuzing 

The  endgame  model  has  the  capability  to  simulate  a conical 
pattern  radar  fuze  with  a back-up  contact  fuze.  The  user  sup- 
plies as  input  to  the  model  tne  dud  probabilities  for  each  of 
these  fuzes.  Thus,  for  example,  by  reading  in  a dud  probabil- 
ity of  1 for  the  contact  fuze,  it  is  possible  to  simulate  a 
missile  with  radar  fuze  only  and  no  back-up. 

2.3.1  Conical  Radar  Fuze 

The  primary  fuzing  mechanism  is  assumed  to  be  a conical  pattern 
radar  device.  It  is  simulated  in  the  endgame  model  by  a right 
circular  cone  whose  apex  is  at  the  missile's  center  of  gravity. 
The  fuze  cone's  axis  coincides  with  the  missile's  body  axis 
and  its  half- angle,  e,  and  range  R,  are  input  quantities.  If 
9 is  read  in  as  exactly  90°,  the  simulation  is  bypassed  and 
fuzing  occurs  at  the  point  of  closest  approach  of  the  missile 
to  the  target. 

In  order  to  analyze  the  simulation  the  missile  body  coordinate 
system  is  convenient.  Thus,  in  missile  body  coordinates  the 
equation  of  the  fuze  cone  is : 

Y2  + Z2  ■=  X2tan2e  . 

The  target  is  represented  by  a set  of  glitter  points.  These 
points  move  relative  to  the  fuze  cone  in  this  coordinate  system. 
Let  the  missile  body  coordinates  of  a particular  glitter  point 
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at  the  end  of  the  previous  interval  be:  Pj : (XG(1) ,YG(1) ,2G(1)) 

and  its  coordinates  now  be:  P2 : (XG(2) , YG(2) , ZG(2) ) . Then,  the 

line  joining  these  two  points  is  represented  by  the  equation: 

XG(1)  - X3  YG(1)  - Y3  ZG(1)  - Z3 

XG( 1)  - XG(2)  " YG( l)  - YG(i)  = ZG(1)  - ZG(2)  = 6 

where  5 locates  any  third  point,  P3 : (X3 ,Y3,Z3)  along  the  line. 

The  point  P3  lies  between  Pi  and  P2  if  and  only  if  0 < 6 < 1. 
Now,  for  P3  to  lie  on  the  fuze  cone,  (X3,Y3,Z3)  must  also 
satisfy : 

Y§  + Z\  - X^  tan2  6 = 0 

Solving  these  two  equations  simultaneously  for  X3  results  in 
the  quadratic: 

AX?  + BX3  + C = 0 

where , 

A = (YG(1)-YG(2))2  + (ZG(1) -ZG(2) ) 2 - tan2 e ■ (XG(1) -XG(2) ) 2 

B = 2 • [(YG(2) -XG(1)-YG(1) -XG(2)) • (YG(1)-YG(2))  + 

(ZG(2)  -XGCD-ZGCl)  • XG ( 2 ) ) ■ (ZG(1) -ZG(2) ) ] 

C = (YG(2) -XG(1)-YG(1) -XG(2)) 2 + (ZG(2) • XG(1) -ZG(1) • XG(2) ) 2 
and  therefore, 

-Bi  Vb2  - 4-A-C 


If  the  roots  are  complex,  the  line  Pi  P2  doesn't  intersect  the 
cone  and  fuzing  does  not  occur  due  to  this  glitter  point  during 
this  time  pulse.  If  the  half-angle  of  the  fuze  is  less  than 
90°  only  positive  roots  are  of  interest.  If  the  half-angle  is 
greater  than  90°  only  negative  roots  are  of  interest.  However, 


in  order  to  simplify  the  computer  coding,  a check  is  made  at 
the  very  beginning  of  the  fuze  subroutine  and,  if  the  half- 
angle  is  greater  than  90°,  the  x- coordinates  of  the  glitter 
points  are  multiplied  by  -1.  Because  the  two  nappes  of  a cone 

% 

are  symmetrical,  this  has  the  effect  of  requiring  all  roots 
of  interest  to  be  positive. 

The  root  resulting  from  the  plus  sign  of  the  radical  represents 
the  earlier  intersection.  This  root  is  calculated  first  and 
tested  to  see  if  it  lies  between  0 and  R-cose.  If  it  is  nega- 
tive then  the  other  root  must  also  be  negative  and  fuzing  does 
not  occur  due  to  this  glitter  point  for  this  time  pulse.  If 
it  is  positive  but  greater  than  R-cose,  then  the  other  root  is 
tested.  If  the  root  does  lie  between  0 and  R-cose,  then  6 is 
calculated  from  the  relationship: 

XG(1)  - X3 
6 = XG(l)  - XG(2)  • 

If  this  quantity  is  greater  than  1 then  it  must  be  greater  than 

1 for  the  other  root  also,  and  fuzing  does  not  occur  due  to 

this  glitter  point  for  this  time  pulse.  If  6 is  negative,  then 

the  other  root  is  examined.  If  0 <_  6 < 1,  then  fuzing  does 

occur  due  to  this  glitter  point  for  this  time  pulse. 
t 

All  glitter  points  are  tested  in  the  above  fashion.  If  none  of 
them  result  in  a value  of  6 between  0 and  1,  then  fuzing  does 
t not  occur  for  the  time  pulse.  Otherwise,  the  minimum  non-nega- 

tive 5 is  used  to  determine  the  time  of  fuzing. 

• 
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2.3.2  Contact  Fuze 

In  order  to  simulate  the  back-up  contact  fuze,  the  aircraft 
is  represented  by  a set  of  direct-hit-without-kill  ellipsoids 
(which  represent  the  skin  of  the  aircraft) . If  the  missile 
trajectory  intersects  any  of  these  ellipsoids  before  the  war- 
head has  detonated,  then  the  contact  fuze  is  initiated.  Thus, 
for  each  of  the  direct-hit-without-kill  ellipsoids  the  follow- 
ing system  of  equiations  in  aircraft  body  coordinates  must  be 
solved . 

Let  the  prior  position  of  the  missile  be:  PI : (XC(1) , YC(1) , ZC(1) ) , 
its  current  position  be:  P2 : (XC (2) , YC (2) , ZC (2) ) , and  a third 
point  along  the  trajectory  be:  P3  : (X3 , Y3 , Z 3)  . If  P3  is  also  on 
the  ellipsoid,  it  must  simultaneously  satisfy: 

XC(1)-X3  YC(1)-Y3  ZC(1)-Z3 

XC(l)-XC(2)  = YC (1) -YC (2)  = ZC(1)-ZC(2)  = d 

and, 

(X;-XO) 2 + (Y3-YO)2  + (Z?-ZO)2  _ 

(XS)2  (YS)2  (ZS)2 

where  XS , YS , and  ZS  are  the  three  semi-axes  of  the  ellipsoid 
and  (XO,YO,ZO)  is  the  center  of  the  ellipsoid. 

Solving  this  system  for  Z3  results  in  the  quadratic,  where: 

A - [(ZC(1)-ZC(2)> -XS-YS]2+t(YC(l)-YC(2)) -XS-ZS]2+ 

[XC(1) -XC(2) ) -YS-ZS]2 

B - -2- {(XC(1)-XC(2)) -YS2 -ZS2- [XC(1) -ZC(2)-XC(2) -ZC(1)+ 

X0- (ZC(1)-ZC(2)) ] + 

(YC ( 1) -YC(2) ) -XS2-ZS2 • [YC(1) -ZC(2)-YC(2) -ZC(1)+ 

Y0- (ZC(1)-ZC(2)) ]+ 

Z0- [XS • YS • (ZC(1) -ZC(2) ) ] 2 } 

::  — —A 
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C = {XS2 -ZS2 • [YC(1) -ZC(2)-YC(2) -ZC(1)+Y0- (ZC(1)-ZC(2) ) ]2  + 

YS2 -ZS2 • [XC(1) • ZC(2) -XC (2) -ZC(1)+X0- (ZC(1) -ZC(2) ) ] 2+ 

XS2 -YS2 • (ZC (1) -ZC (2) ) 2 • (Z02-ZS2) } 

and  therefore, 

„ _ -B±  Vb2  - 4-A-C 

2 a 

If  the  roots  are  complex,  no  direct  hit  occurs  on  this  ellip- 
soid for  this  time  pulse.  If  real  roots  do  exist,  the  root 
which  represents  the  earlier  intersection  is  the  root  arising 
from  the  case  where  the  radical  has  the  same  sign  as  (ZC(1)- 
ZC(2)).  This  root,  therefore,  is  the  only  one  that  is  examined. 
The  quantity  6 is  calculated  from  the  relationship: 

ZC(1)-Z3 
: = ZC ( 1) -ZC (2) 

If  0 < { < 1,  then  P3  lies  between  Pj  and  P2 , and  the  solution 
represents  a valid  direct  hit  for  the  current  time  pulse. 

All  direct-hit-without-kill  ellipsoids  are  tested  in  this  manner 
and  the  solution  which  gives  rise  to  the  least  positive  6 is 
used  to  calculate  the  time  at  which  the  contact  fuze  is  initiated. 

2 . 4 Single-Shot  Component  Kill 

Three  kill  mechanisms  are  considered  in  the  endgame  model. 

These  are  direct-hit  kill,  blast  kill  and  fragment  kill.  Direct- 
hit  kill  occurs  whenever  the  as  yet  undetonated  missile  strikes 
the  aircraft  in  such  a way  as  to  cause  a kill  even  without  the 
warhead  effects.  Blast  kill  results  from  the  shock  wave  formed 
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in  the  surrounding  atmosphere  upon  detonation  of  the  explosive 
warhead.  This  shock  wave  can  damage  an  aircraft  at  close  range 
by  overpressure  and  impulse  loading.  Fragment  kill  occurs  when 
fragments  of  the  warhead  casing,  propelled  by  the  force  of  the 
warhead  detonation,  intercept  one  or  more  of  the  vulnerable 
components  of  the  aircraft. 

2.4.1  Direct-Hit  Kill 

Direct-hit  kill  is  simulated  in  the  endgame  model  by  represent- 
ing the  target  as  a set  of  direct-hit-with-kill  ellipsoids.  The 
process  of  determining  whether  or  not  such  a direct  hit  occurs 
in  particular  time  pulse  is  identical  to  the  process  described 
in  Section  2.3.2  except  that  the  direct-hit-with-kill  ellipsoids 
are  used  in  place  of  the  direct-hit-without-kill  ellipsoids. 

2.4.2  Blast  Kill 

A balst  kill  occurs  whenever  the  burst  point  (X*,Y*,Z,V)  lies 

within  one  or  more  of  a set  of  blast  ellipsoids.  The  equation 

expressing  this  condition  in  aircraft  body  coordinates  is : 

(X*-XBL) 2 + (Y*-YBL) 2 + (Z*-ZBL) 2 < l 
(ABL)2  (BBL) 2 (CBL) 2 

where  (XBL, YBL.ZBL)  is  the  center  point  of  a blast  ellipsoid  and 
ABL,  BBL  and  CBL  are  the  scaled  semi-axes  lengths  for  that  ellip- 
soid. This  test  is  performed  for  all  blast  ellipsoids  and  blast 
kill  is  said  to  occur  if  it  is  satisfied  for  any  one  of  them. 

2.4.3  Fragment  Kill 

The  static  explosion  of  a fragment  warhead  produces  a character- 
istic spectrum  of  fragment  mass,  fragment  density  and  fragment 


emission  speed.  The  explosion  of  a miving  fragment  warhead 
alters  this  spectrum  by  virtue  of  the  imposed  forward  speed 
of  the  missile.  It  is  necessary  to  determine  the  interaction 
of  this  altered  spectrum  with  the  moving  aircraft  in  order  to 
obtain  expressions  for  the  expected  number  of  fragment  hits 
upon  the  vulnerable  components  of  the  aircraft  and  for  the  net 
striking  speed  of  these  fragments. 

The  distance- time  equation  governing  the  fragments  can  be 
written : 

ln(l+CD-m"1/3-o-V0-t)  = CD-m-1/3-P-L 

where,  is  the  sea  level  drag  coefficient  for  the  fragment, 
m is  the  fragment  mass,  p is  the  relative  air  density,  V0  is 
the  initial  fragment  speed,  t is  the  time  of  flight  of  the 
fragment  and  L is  the  distance  travelled  by  the  fragment.  This 
equation  has  no  closed  form  solution.  The  endgame  program 
arrives  at  a solution  using  the  Newton-Raphson  iterative  method 
to  find  a root  of  the  equation: 

F(t)  = ln(l+CD-m"1/3-p -Vo-t)  - CD-m‘1/3-p-L  = 0 

This  classical  method  takes  the  Jth  estimate  of  the  root,  t^3\ 
and  extracts  the  (J+l)th  estimate,  t^3+^^ , by  means  of  the 
following  equation: 

r(J+l)  _ rU>  F(t(J)) 
t = t - 

F(tU;) 

The  procedure  is  repeated  until  successive  estimates  are  con- 
sidered to  differ  neglibibly. 
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As  a first  estimate  of  the  solution  is  used  the  time  of  flight 
resulting  from  the  case  where  the  fragment  does  not  slow  down 
due  to  drag.  The  solution  for  the  fragment  time  of  flight  and 
distance  travelled  also  determines  the  dynamic  emission  angle, 

Y » between  the  missile  velocity  vector,  V^,  and  the  dynamic 
fragment  velocity  vector,  Vq.  This  in  turn  determines  the  angle 
e,  between  the  static  fragment  velocity  vector,  V£ , and  (see 
Figure  2-2) . 

As  a result  of  the  rotation  of  velocity  vectors  due  to  the  mis- 
sile motion,  the  fragment  density  in  the  static  case,  H' , is 
unequal  to  the  dynamic  density, 

’dyh  - ,/E  ■ 

The  term  E can  be  shown  to  be : 

E = s^-n  y . d'l 

sine  3e 

A2 • jA+cose ] 

[A::+2  • A • cosS+1  ] 3/2 


where , 


VVM 


Input  data  exist  for  the  fragment  density  resulting  from  a 
static  explosion  in  the  form  of  a table  having  entries  at  10° 
intervals  off  the  nose  of  the  missile.  Thus,  in  order  to  in- 
terpolate for  ij-  from  the  input  table,  it  is  necessary  to  find 
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into  missile  body  coordinates  using  the  matrix  TAM  and  taking 
the  dot  product  of  this  with  the  vector  (1,0,0). 

In  the  Dynamic  SAM- Aircraft  Model  the  aircraft  velocity  vector, 
VT>  is  not  necessarily  parallel  to  the  aircraft  body  axis. 

Thus,  in  the  aircraft  body  coordinate  system: 

VT  = VT(l)-i  + Vt(2)*J  + VT(3)-k 

Consider  now  a coordinate  system  which  is  identical  to  the  air- 
craft body  system  at  the  time  of  detonation  but  which  remains 
fixed  in  space.  Let  (X0,Y0,Z0)  be  the  coordinates  of  one  of 
the  vulnerable  components  of  the  aircraft  in  this  stationary 
system  at  the  time  of  detonation.  Then,  if  TF  is  the  time  of 
flight  of  the  fragment,  the  coordinates  of  this  component  at 
the  time  it  is  hit  will  be: 

Xjj  = X0  + VT(1)  TF 

YH  = Y°  + V2)  TF 
ZH  = zc  + Vt(3)  TF 

As  previously  defined  in  Section  2.4.2,  (X*,Y*,Z*)  represents 
the  coordinates  of  the  burst  point  in  the  aircraft  body  system 
at  the  time  of  detonation.  Therefore, 

L - VCXy-X*)*  + (Yh-Y*)"'+  <Zh-Z*)^ 


and  the  direction  cosines,  with  respect  to  the  aircraft  body 
system  axes,  of  the  line  directed  from  the  explosion  point  to 
the  hit  point  are  Sjj,  6y  and  e^: 


r 
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cT-x* 


"X 


Y -Y* 
H 


bY  = L 


vz* 


"Z 


The  above  method  of  calculating  hits  assumes  that  there  are 
no  masks  or  shields  to  stop  the  fragments  before  they  reach 
the  vulnerable  component . Shielding  is  incorporated  into  the 
model  by  the  approximate  but  adequate  treatment  described 
below. 

The  hit-point  location  of  a vulnerable  component  is  calculated, 
assuming  that  the  shield  is  non-existent.  This  defines  the 
time  of  flight  of  the  fragment,  TF,  and  also  the  straight  line 
path  of  the  fragment  (from  explosion  point  to  hit  point) . The 
shield  is  now  placed  in  the  position  it  would  occupy  at  the 
time  of  the  hit.  The  shield  is  taken  to  be  an  ellipsoid  whose 
axes  are  parallel  to  the  X,  Y,  Z axes  of  the  aircraft  body  co- 
ordinate system.  If  the  fragment  travel  line  intersects  the 
shield,  the  vulnerable  component  is  considered  to  be  shielded. 
If  there  is  no  intersection,  the  component  is  not  shielded. 

This  treatment  should  be  accurate  if,  at  the  time  of  the 
explosion,  the  shield  is  much  closer  to  the  component  than  to 
the  explosion  point.  A higher  level  of  accuracy  would  involve 
ballistic  solutions  for  fragment  and  shield.  This  appears  both 
complicated  and  unwarranted. 
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TABLE  2,1 

CRITERIA  FOR  ASPECTS  STRUCK 


Criterion 


Aspect  Struck 


VT(1) 

Sy  > y 

A VHIT 


Rear 


VT(1) 

ex  " \T 

X VHIT 


Front 


Vt(2) 

• y > V 
1 VHIT 


Right 


Y < V, 


Vt(2) 


Left 


Vt(3) 
°Z  > V 

^ 'hit 


Bottom 


Vt(3) 
< V 

L HIT 


To  calculate  the  number  of  fragment  hits  on  the  component, 
it  is  necessary  to  know  the  fragment  density  along  the  vector 
VNET'  aPPears  to  be  very  difficult  mathematically  to 

arrive  at  this  density,  however,  Therefore,  the  approximation 
is  made  that  the  number  of  hits  can  be  calculated  on  a static 
target  (equivalent  to  using  the  density  ^Dyn'  a^on6  t*ie  frag- 
ment  velocity  vector).  This  is  an  excellent  approximation  if 
Vrit  >>  V^,.  In  situations  in  which  this  condition  does  not 
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hold,  the  striking  energy  of  a fragment  would  be  so  low  as 
not  to  cause  an  appreciable  probability  of  kill. 


Using  the  static  target  concept,  then,  Table  2.1  is  modified, 
replacing  the  terms  VT(1)/VHIT,  VT(2)/VHIT  and  VtO)/vHiT  by 
zero.  The  vulnerable  areas  of  the  relevant  aspects  are  pro- 
jected onto  a plane  perpendicular  to  the  fragment  velocity 
vector  in  the  following  manner: 


VA  = 


•VA 


X 


>Y 


•VAy 


■VA, 


The  expected  number  of  lethal  hits  on  the  component,  N,  can 
now  be  calculated: 


N 


DYN 


•VA/L2 


Applying  Poisson's  Law  gives  an  expression  for  the  probability 
of  survival  of  the  component: 


The  above  calculations  are  carried  out  in  the  endgame  program 
for  all  components  and  at  all  kill  levels. 

2 . 5 Overall  Aircraft  Kill 

The  single-shot  aircraft  kill  probability,  PKAC,  is  calculated 
from  the  following  relationship  for  an  aircraft  composed  only 
of  singly  vulnerable  components. 

NCOMPS 

PKAC  - 1 - n ps. 

i-1  1 


' » . - ....  . 
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where  PS^  is  the  survival  probability  of  component  i and 
NCOMPS  is  the  total  number  of  vulnerable  components. 

This  formulation  is  altered  somewhat  when  there  are  one  or 

more  pairs  of  doubly  vulnerable  components.  If  j and  k form 

a pair  of  doubly  vulnerable  components,  then  the  term  PS. -PS, 

3 k 

is  replaced  in  the  above  product  by  the  probability  that  both 

j and  k are  not  killed:  [ 1- (1-PS  . ) • (1-PS,  ) ] . 

J K 

In  order  to  calculate  the  cumulative  aircraft  kill  probability 
resulting  from  a salvo  of  missiles  it  is  necessary  to  compute 
and  save  the  cumulative  survival  probability,  CPS^,  for  each 
vulnerable  component,  i : 


cps.  = n PSi 

all 

shots 

The  cumulative  aircraft  kill  then  becomes 


CPKAC 


NCOMPS 

n 


cps1 


with  CPSj-CPSk  being  replaced  by  [1- (1-CPSj) • (1-CPS^) ] for 
any  pair,  j and  k,  of  doubly  vulnerable  components. 
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